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1. Introduction 

In Paper I fll], we derived exact properties of the Einstein-Boltzmann equations with 
a relaxation-time model of collisions. An exact truncated distribution was used to 
derive transport equations similar in form to the Israel-Stewart thermodynamics. We 
also found an exact truncated solution for massless particles in a flat FRW spacetime. 
This paper extends Paper I by considering exact non-truncated Einstein-Boltzmann 
solutions in FRW and Bianchi I spacetimes. 

In Section 2 we present exact anisotropic solutions for a flat FRW spacetime, and 
investigate their properties. In particular we show that the bulk viscosity of these 
solutions is always zero, since they are relaxing towards free-streaming equilibrium. 
For solutions that relax toward hydrodynamic behaviour, we give an exact formula for 
the bulk viscous pressure. Our solutions generalise the equilibrium anisotropic solutions 
given by . Isotropic non-equilibrium solutions are only possible for massless particles, 
and they have the surprising property that entropy is generated despite the vanishing 
of bulk viscosity, heat flow, shear viscosity - and of all non-scalar moments of the 
distribution. This illustrates the point that the standard dissipative quantities (bulk 
viscosity, heat flux, shear viscosity) cannot provide a complete or exact description of 
non-equilibrium states. 

In Section 3 we find exact anisotropic solutions in Bianchi I spacetimes. The 
solutions have zero heat flux and bulk viscosity but nonzero shear viscosity. In Section 
4 we use the results of the previous sections to construct a model in which an Einstein- 
Boltzmann solution in FRW spacetime evolves into an Einstein-Boltzmann solution in 
Bianchi spacetime. The basic idea, due to Matravers and Ellis 0, is that anisotropy in 
the FRW distribution is communicated to the geometry via the Einstein field equations 
when the gas cools sufficiently for massive particles to become non-relativistic. We 
provide an explicit dissipative realisation of the Ellis-Matravers model, which is in 
accord with Penrose's Weyl Curvature Hypothesis, i.e. that the universe is initially 
conformally flat and that anisotropy, inhomogeneity and entropy production develop as 
the universe expands 0,0. 

The distribution f(x,p) satisfies the Boltzmann equation with relaxational collision 
term, i.e. the BGK equation ]I[ 

where 7 encodes microscopic interaction information in a linear, macroscopic 
approximation |l[], / is the distribution towards which (or away from which if 7 < 0) / 
is relaxing, and E = —Uip 1 is the particle energy relative to the average four-velocity u l 
associated with /. If / is a dynamic (or 'global') equilibrium distribution |TJ, L[f] = 0, 
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then the solution of flip is 

/= f + he~ r where L[h] = 0, T (x(v), E(v)) = J* 7 (x(u), E{u)) du (2) 
The standard form for 7 is the Anderson- Witting (AW) form 0: 

where r is the mean interaction time. This form includes the case of radiative transfer 
with isotropic scattering Using dv = dt/p° = dt/E and ([3p, the relaxation factor in 
(0) becomes a spacetime scalar: 

The distribution may be covariantly decomposed relative to « ! M: 

f(x,p) = F(x,E) + F i (x,E)e i + F ij (x,E)e i e j + ■■■ (5) 

where the covariant multipoles F^... are isotropic, spatial, trace-free and symmetric, and 
e l is the unit spatial projection of p l . 

In order to avoid unnecessary details, we will not summarise the basic equations 
and results of Paper I, but refer where necessary directly to the equations in that paper 
in the form (I:n), where n is the equation number in Paper I. 

2. Exact non— equilibrium Einstein— BGK solution in FRW spacetime 

In Paper I we found an Einstein-BGK solution for a truncated form of the distribution 
function. This solution was used to derive a set of exact thermodynamic laws. In this 
section we present the full non-truncated solution in flat FRW spacetime with natural 
coordinates x % = (t,x v ): 

ds 2 = -dt 2 + R 2 (t)(dx 2 + dy 2 + dz 2 ) (6) 

By homogeneity, the spatial momenta p v are constants of the motion, i.e. L[p u ] = 
0. Thus homogeneous distributions are of the form f(t,p u ), while isotropic and 
homogeneous distributions are of the form f(t,w) where 

w 2 = ( Pl ) 2 + (p 2 ) 2 + (p 3 ) 2 = (E 2 - m 2 )R 2 (7) 

and where u l = 5 l $. If the homogeneous / is relaxing toward the homogeneous and 
isotropic equilibrium, /, then by (§) and @ 

ft di 1 

f(t, Pu ) = f(w) + %,)e- r « , r(f) = / — - (8) 



/o r(t') 

Clearly / is spatially homogeneous but dynamically anisotropic, and it must depend 
explicitly on the cosmic time (i.e. T 7^ 0) if it is non-equilibrium. Note that since 
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L[f] = 0, / can only be a collision-dominated Maxwell-Boltzmann equilibrium if 
m = 0; for m > 0, / must be a collision-free equilibrium distribution [DJ. In this case 
the BGK solution / represents a distribution relaxing towards free-streaming isotropic 
equilibrium. 

In order to investigate the properties of the solutions, we use the covariant harmonic 
decomposition (^: 

F{t, w) = F(w) + H(w)e- r ^ , F u ... K (t, w) = H v ... K {w)e- T ® (9) 

The kinematics and the dynamics of the solution (§) are determined by the particle 
four-current N % = nu l + k % and energy-momentum tensor = fiitiitj + phij + 7Tjj + 
qiUj + qjUi, where n is the number density, k l is the number flux, \x is the energy density, 
p is the isotropic pressure, = + u(Uj is the spatial projector, ir^ is the anisotropic 
pressure tensor and qi is the heat flux. These quantities follow from (1:40) and (^): 

n =wH / u> H(w)dw (10) 

-R (t) Jo 

^ = ^77r/ ^K + m 2 ^ 2 ^)]- 1 ^^)^ (11) 
SK^it) Jo 

£ + -J^TT" / W V + m 2 ^ 2 ^)] 1 / 2 ^^)^ (12) 
it (rj jo 

P+^WI wV + rn 2 R 2 (tT 1/2 H(™)dw (13) 

3i? 4 (i) 
87re- r « 



P 



7T„ 



w A H v (w)dw (14) 

/•CO 

l r )/?'(/) /o w V + m W^H4» (is) 

with k = 0, go = and 7r j = 0. A non-zero number flux gives a particle drift that 
is out of keeping with FRW symmetry, and although it is possible to satisfy the field 
equations for fc 8 ^ (see the fluid solutions of Calvao and Salim 0), we regard this 
as unnatural. From equation (11) it is clear that in order to get a zero number flux, 
which gives a non-tilted kinematic average 4-velocity, H v must vanish for m > 0. It is 
possible to find nonzero H v if m = 0: 

/ w 2 H u (w)dw = (16) 
Jo 

The full Boltzmann collision term is based on microscopic conservation so that the 
macroscopic conservation of momentum and energy are identically satisfied. This is not 
the case for the BGK collision model, and the conditions imposed by the conservation 
equations require separate investigation QTJ| . The conservation of particle number, energy 
and momentum are given by equation (1:48). On using the FRW BGK solution (P), we 



5 



find the following condition for the conservation of particles: 

roo 

/ w 2 H{w)dw = (17) 
Jo 

which by (10) implies n = n, i.e. the number density is matched to that of the limiting 
equilibrium. 

The condition for the conservation of energy gives 

w 2 [w 2 + m 2 R 2 (t)} 1/2 H(w)dw = (18) 

If m > 0, this forces H = 0, but H may be nonzero for m = 0. In all cases, ( |T8| ) in (12) 
implies \i — p,, so that the energy density matches the equilibrium value. 
Momentum conservation gives 

POO 

/ w 3 H u {w)dw = (19) 
Jo 

Condition fllTf ) is automatically satisfied for m > 0, since H u = from k u = 0. For 
m = 0, it is a further condition. By (14), we see that (|i~9"|) leads to g y = = 0, so that 
is also the energy-frame four-velocity (it is already the particle-frame four-velocity 
by k v = 0, so the two four-velocities coincide in these solutions). 

The entropy density s is given by (I:40g). For the FRW solution (|9|): 



47T 



R 3 (t) Jo 



roo r r — \i r — 

J w 2 { [F(w) + H^e^l [l - ln[F(w) + H(w)e 



-r(*)i 



F(w) + H(w)e 



-r(t) 



e- 2m HJw)H u (w) 



dw 



(20) 



The entropy production rate (I:40i) becomes 

° {^(w) In [FH + H{w)e- r{t) 



S' 



w 



R 3 (t)r(t) Jo 

T{t) H u (w)H v (w 
2F(w) + H(w)e- m 



x 



F(w) + H(w)e~ m 
+ • ■ - j dw; 



x 



(21) 



For the Boltzmann collision term, S l -i > follows identically, but this is not true for 
the BGK collision term [|TJ. It is not clear in general from (21) whether the H-theorem, 
S l ; i > 0, is satisfied without restrictions for the BGK solution (§). This needs to be 
checked for each specific solution. Note that for m > 0, when H = = H u , the only 
contribution to S" 1 ^ is from the quadrupole and higher moments. For m = 0, if is in 
general nonzero and there is a monopole contribution to the entropy production: 

47re- r « 



In summary: 

H = 



R 3 (t)r(t) Jo 



w 2 H(w) \n[F(w) + H(w)e~ l w ]dw + 



for m > 



roo roo 

/ w r H(w)dw = 0= / w r H u (w)dw for m = (r = 2,3) 
Jo Jo 



(22) 

(23) 
(24) 
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With (23) and (24), conservation of particle number and energy-momentum is satisfied 
and we have 



n 



n, n = ji , k u = , q v 







(25) 



Now we impose the Einstein field equations (1:61) (note that the conservation equations 
are already satisfied): 



= lt VK 

3^ 
6 R? 



(26) 
(27) 



The heat flux is already zero by (^5|) . Vanishing anisotropic stress requires, by (15) 



4r 2 
W \W 



+ m 2 R 2 (t)}- 1/2 H VK (w)dw = 
For m > (|28| ) forces FL VK = 0. For m = 0, H UK is subject to 

POO 

/ w 3 H UK (w)dw = 
Jo 

Using (12) in the Friedmann equation (27) we can write it as 

Air 



(28) 



(29) 



R(t) 



/ w 2 [w 2 + m 2 R 2 {t)] 1/2 F(w)dw 
Jo 



3R 2 (t) 



47re" r(i) r°° 



n l/2 



By (23) and (24) the term containing H(w) in (30) vanishes for m > and we can 
re-arrange this equation and give the solution explicitly as 



V3 f R2 



Air 



w 2 [w 2 + m 2 u] 1 / 2 F(w)dw 



-i/2 



du 



(31) 



which is the same as for an equilibrium Einstein solution J|. By specifying F, R(t) can 
be determined in principle, and the metric @ will be known - completing the Einstein 
solution. This is analogous to specifying an equation of state in a fluid model. 

The reason that R has the same form as for an equilibrium solution lies in the 
vanishing of the bulk viscous pressure U = p — p. By (13) 

Attp~ r r°° 

{62) 



n 



3i? 4 



4 [ W 2 + m 2 J R 2 ]~ 1 / 2 ^Hrf 



w 



It follows from (23) and (24) that II = for m > 0. For m = 0, this is in accord 
with the approximation schemes used to derive transport equations [0. For m > 0, 
the vanishing of bulk viscosity is a consequence of our choice of / satisfying L[f] = 0. 
This means that / cannot be a Maxwell-Boltzmann distribution (unless the universe 
is static) P[, and therefore / cannot be approaching the hydrodynamic regime where 
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the standard approximation schemes are applied and predict II 7^ 0. Instead / must be 
collision-free, and if / is near to free-streaming , our results show that, unsurprisingly, 
there is no bulk viscosity. 

If we drop the restriction L[f] = and take / to be a local Maxwell-Boltzmann 
distribution, we have M 

E (33) 



/ = exp 



a(t) 



T(t)_ 

where T is the temperature and a the chemical potential. Using the general (i.e. 
L[f] 7^ 0) BGK solution (1:44), we find the following exact formula for IT: 



IT = - — e" r / " dt e r+a I " dE (E 2 - m 2 f 2 e- E ' T x 



3 

x I a + 



f 


dt 


e r +a 


Jo 






E 


f 


R 
+ R 


f 


f 



2 ' 

m 



(34) 



In the case m = 0, when T cx and a = 0, we have II = 0, as expected. 

What is remarkable about the solutions with L[f] = is this: despite the vanishing 
of viscosity and heat flux, there is dissipation, since L[f] 7^ 0. The simplest case is the 
isotropic massless solution (a constant) 

f = f + He- T ®, f = e a ~ w (35) 

first presented in ffl]. By fl22|) and (24) the entropy production rate is 



S\i = R , {t)T{t) I w 2 H(w) ln[l + H(w)e^- a ]dw (36) 

and this is clearly positive since / is. The point is that the bulk viscosity, heat flux 
and shear viscosity are only approximate indicators of non-equilibrium states, and 
even if they vanish there can still be dissipation. Non-equilibrium states in general 
cannot be completely, and certainly not exactly, described by these standard dissipative 
quantities. Usually this is understood in terms of the effect of multipoles higher than 
the quadrupole, which are neglected in the standard approximation schemes. However, 



the isotropic solution ( p5|) has no multipoles beyond the scalar monopole, and yet it is 
out of equilibrium. 

3. Exact non— equilibrium Einstein— BGK solutions in Bianchi I spacetime 

Bianchi I spacetime, 

ds 2 = -dt 2 + X 2 (t)dx 2 + Y 2 (t)dy 2 + Z 2 (t)dz 2 (37) 

is distinguished kinematically from FRW spacetime by non-zero shear, whose evolution 
is given by 

<jij - Pike j + \hijV kl &ki - §7!"ij + Eij = (38) 



where E ik = Cijkiu^u 1 is the electrical part of the Weyl tensor Cyjy. By symmetry, the 
trace-free spatial tensors have the form 

A*j = diag (0, A\, A\, -A\ - A 2 2 ) (39) 

where A^ = a^, 7Ty, £y or i^-. 

The Einstein field equations for this metric are || 

-|(/i + 3p) (40) 

(41) 

K/z-pHtt 1 ! (42) 

(43) 



X 


Y 
+ Y + 


z 


X 


z 


Qj 






X 


YX 
+ YX 


zx 
+ ~zx 




Y 


YX 
+ YX 


ZY 
+ ~ZY 


Y 


Z 


ZX 


ZY 


Z 


+ zx 


+ ~ZY 



r 2 2 



1 ^-p)-n 1 1 -n 2 2 (44) 

The conservation of energy-momentum reduces to (0 = u l -i) 

/i + (/i+p)e + 7Ty(7 ij =0 

which is identically satisfied if (40) - (44) are satisfied. 

By homogeneity, the spatial momenta p u are constants of the motion. Thus 
homogeneous distributions are of the form f(t,p u ), including the special case f(t,w) 
where 

w 2 = (pi) 2 + (p 2 ) 2 + (p 3 ) 2 

= X\t){p 1 ) 2 + Y\t){p 2 ) 2 + Z\t)(p 3 ) 2 (45) 

Contrary to the FRW case, w is not isotropic, i.e. it is not a function only of E in 
momentum space, where 

E = [m 2 + X- 2 {t){ Pl ) 2 + F- 2 (t)(p 2 ) 2 + Z- 2 {t){pz) 2 ] 1 ' 2 (46) 

Both w and E are spatially homogeneous. However, there is no simple relation between 
the anisotropic constant of motion w and the isotropic non-constant energy E, unlike 
the FRW case. This happens since hij is anisotropic: 

A 2 = E 2 — m 2 = h^p^p u 
w 2 = 8^ PflPv 

In the FRW case = Rr 2 {t)5» v . 

If the homogeneous, anisotropic / is relaxing toward the homogeneous, anisotropic 
/, where L[f] = 0, then the BGK solution has the same form (^) as in the FRW case, but 
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with w given by (45) (and anisotropic). It follows that the Bianchi I solution matches 
the FRW solution (R|) in the limit R(t) = X(t) = Y(t) = Z(t). The covariant harmonic 
decomposition (|5J) takes the form 

/(*, Pv ) = F(t, E) + F K (t, E)e K + F Kp (t, E)e K e? + ■■■ (47) 

where 

F„... K (t, E) = F v ... K (t, E) + e- m H u ... K (t, E) 

Note that since / is anisotropic, its higher order multipoles cannot be neglected in the 
decomposition. Although the solution and its decomposition do not take the convenient 
form of the FRW case, it gives us the necessary tools to investigate the conditions for 
an Einstein-BGK solution. 

Because of the nonvanishing shear in the Bianchi I geometry, the covariant 
multipoles in fl47|) are no longer independent. The relationship between the multipoles 
is determined by the Boltzmann equation. By attaching an orthonormal tetrad to u l , 
the Liouville operator L can be written in the covariant harmonic form (1:10). This 
allows one to write the Boltzmann equation as a set of coupled differential equations 
in the multipoles. The first two Boltzmann multipole equations for a homogeneous 
distribution function [11] become, for a Bianchi I geometry and a BGK AW collision 
term: 

BF 

-E— = r- 1 E(F-F) (48) 

+ ^A 2 e|| = r-i E (F u - F v ) (49) 

where we have used the Bianchi I symmetry to simplify the expressions given in 0. 
The higher order multipole equations O] (p492) together with (48), (49) show that if 
F, F u are specified, then the multipole equations place a chain of restrictions on the 
quadrupole and higher moments. The important point is that if the shear is non-zero, 
the multipoles are no longer independent. 

From (1:40), the kinematic and dynamic quantities of the solution fl47D are 

n = n 

J m 

- Att r°° 

k v =k u + —e- r{t) \ 2 H u (t,E)dE (51) 



+ 4vre- r(i) / EXH(t, E)dE (50) 

J m 



3 
An 



poo 

--fi + 4ne- m E 2 \H(t, E)dE (52) 

Jm 

Air r°° 

p =p+ —e~ m / X 3 H{t,E)dE (53) 

3 Jm 
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Air 



Air r°° 

Qv =qu + —e~ m E\ 2 H l/ (t,E)dE (54) 

O Jm 

Sty r°° 

Kvn = n VK + —e~ m / X 3 H UK (t, E)dE (55) 

ID J m 

The particle flux ki, energy flux and anisotropic stress 7fjj of the equilibrium 
distribution / are in general nonzero because / is anisotropic. These quantities do not 
reflect any dissipation, but are part of the measure of deviation from isotropy (compare 



0, 0)- 

Given that the BGK collision model does not guarantee macroscopic conservation, 
we impose the conditions for the conservation of particle number, energy and momentum 
(I:40g): 



III 

oo 



EXHdE = (56) 

E 2 XHdE =0 (57) 

E\ 2 H u dE = (58) 

With (50), (52) and (54) these give the matching conditions: 

n = n, fi = ft, q v = q v (59) 

The Bianchi I solution has equilibrium particle number density, energy density and 
energy flux. In general, since H and H u depend explicitly on time, (56) - (58) require 

H = = H„ (60) 

The Einstein field equations are imposed next. Equation (41) along with (54) and 
(58) implies 

POO _ 

E 2 X 2 F u dE = 



L 



Again, in general this requires 

F v = (61) 
Equations (p9|) and (55) require 

/ \ 3 F\{t,E)dE=— diagfTrVTr^-A-TrV) (62) 

Thus Fij has at most two independent components. Then (48), (49), (|59|), (|61|) and 
(|62|) are the restrictions on the harmonics for an Einstein-BGK solution. Once F and 
F v v (no sum) are specified, the remaining field equations in principle yield a solution 
{X(t),Y(t), Z(t)} (this is analogous to specifying equations of state for p and tt^- in 
fluid models). 
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By (53) the bulk viscous pressure IT is given by 
Ait f 00 

II = — e~ m / X 3 H{t, E)dE (63) 

O J m 



It follows from (|59|) and (|60| ) that II = for m > 0. Again the bulk viscosity vanishes 
when m > because / is relaxing toward free-streaming. 



The restrictions (|60|), (|61| ) and ( |62D on the multipoles may for example be satisfied 
by the choice: 

F = F, F u = = F u , F WK = 0, H\ = dmg(V 1 ,V 2 ,-V 1 -V 2 ) (64) 

where 

ir\ = —e- m \ 3 V„(t,E)dE (no sum) 

15 Jm 

This choice implies that the equilibrium distribution / has perfect fluid behaviour, i.e. 



K = 0, q v = , tt„ k = 

and by (51), there is no particle flux, i.e. k u = 0. 

The BGK multipole equations (48) and (49) reduce to 

2e ~ T jE ( AV ^-) + 5A30 § - 15XE % = ( 65 ) 
<} {\ 4 a Kp F UKp ) =0 (66) 



BE 

The higher order BGK multipole equations become conditions on the fourth and 
higher multipoles which may always be satisfied, since these multipoles are otherwise 
unrestricted. The condition (66) gives 

o Kp F VKp = 

while (65) is a constraint on F and V v . Once these are specified subject to (65), the 
remaining field equations determine g vv {t) in principle, thus completing the Einstein- 
BGK solution. Using the fact that 

^W§_ie|-ie,f-ie) 



(65) leads to the new result 

- r (!-I)^(^) + - r (M)^(^ 

It is clearly possible to find F, V v that satisfy the single linear equation (67). 
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4. Anisotropy generation in FRW cosmologies 

Using the FRW and Bianchi I results we now construct a model in which an Einstein- 
Boltzmann solution in FRW spacetime evolves into an Einstein-Boltzmann solution in 
Bianchi I spacetime. It is assumed that the universe initially has FRW geometry with 
a matter distribution that is described by an anisotropic distribution function which 
is also compatible with FRW geometry. Ellis and Matravers fl4] , Hl5| , || propose two 



mechanisms whereby the anisotropy of the matter distribution is communicated to the 
universe. The first applies to a universe that is initially FRW with an equilibrium particle 
distribution that is inhomogeneous and anisotropic. The anisotropy and inhomogeneity 
is not communicated to the spacetime geometry. The particle distribution is effectively 
collision-free, because the particles are assumed to enjoy asymptotic freedom under 
the high temperature conditions of the early universe. As the universe expands, the 
temperature drops and the collisions become significant, allowing the anisotropy to be 
communicated to the geometry. In the second mechanism, while the temperature is 
very high the particles effectively behave as if their rest mass is zero. As the universe 
expands and cools, and the threshold energies of different particles are reached, their 
rest mass becomes significant. Again, this effective change in the zero rest mass allows 
the inhomogeneity and anisotropy to be communicated. We show that the generation 
of anisotropy can be modelled starting from an anisotropic and non-equilibrium BGK 
solution. A non-equilibrium model allows us to overcome some of the drawbacks of 
the equilibrium model - for example, the difficulty in motivating a collision-free early- 
universe phase. 

We assume that the universe initially (t < to) has FRW geometry and that the 
particle distribution is given by (H). This distribution function is dynamically anisotropic 
and has Bianchi I symmetry, i.e / = f(t,p v ). For the high temperature conditions in 
the early universe, the particle rest mass is insignificant and an Einstein-BGK solution 
may be chosen that is consistent with the Bianchi I solution fl64l) . Furthermore, the 
local equilibrium / to which / is close is a collision-dominated Planckian distribution, 
in keeping with standard physics of the early universe (and not requiring asymptotic 
freedom or other exotic processes). 

By (pD, (24), (H) we can take 

r OO 

H = 0, H v = 0, / w 3 H UK {w)dw = (68) 
Jo 

where w = R(t)E. Thus under high temperature conditions for which m — 0, it is 
possible to find distribution functions with non-zero quadrupole, but for which the 
condition ir UK = is satisfied. We need non-zero H VK in order to 'switch on' n VK (the 
mechanism for this will be described below) . Then by fl38| ) the shear anisotropy <Jij will 
emerge and the geometry will evolve to a Bianchi I phase. By ( p>3| ) and (55), we could 
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take H VK of the form 

H" K (w) = diag [Vi(w), V 2 (w), -V^w) - V 2 (w)} (69) 
where V v are non-zero but satisfy 

w z V u (w)dw = Q (70) 

o 

Thus for t < t Q the universe is considered FRW, while for t > to it is Bianchi I. On the 
hypersurface t = t Q , the geometry and the distribution function symmetry must satisfy 
both the conditions for FRW and Bianchi I spacetimes. This is the case if the following 
matching conditions are satisfied 0: 

X(t ) = Y(t ) = Z(t ) and X(t ) = Y(t ) = Z{t ) (71) 

As discussed in 0] and [|15||, these matching conditions always have a solution. Then 



this solution becomes the initial conditions for the Bianchi I solution, which is governed 
by existence and uniqueness theorems for the Einstein-Boltzmann equations (see ||15||). 
The transition from effectively massless to effectively massive behaviour will take place 
over a cosmic time 5t which is much less than the expansion time. Evolution away 
from FRW is triggered as soon as the mass becomes dynamically significant, i.e. after 
time St. In practice we are treating the transition as instantaneous, in a similar way to 
standard models of electron-positron annihilation. Our simple model illustrates that in 
principle, finely tuned anisotropy in the matter distribution can be unlocked dynamically 
to generate anisotropy in the spacetime geometry. 

Note that the evolution away from FRW occurs entirely within the Bianchi I phase, 
as the shear grows from zero. At the transition, 0"jj(to) = 0, and <Jij remains zero unless 
TTij becomes non-zero to force the universe to evolve to the Bianchi I geometry. In 
the Bianchi I phase t > t , the solution is given by (|47j ) and it clearly reduces to the 
FRW solution at t = t as a result of the matching conditions ([71]) . Because w becomes 
isotropic and reduces to the FRW form with the application of the matching conditions, 
the Bianchi I equilibrium solution / reduces to the FRW equilibrium solution at t . 

In summary, the mathematical model is the following. The collision model is AW: 

rt dt' 

T(t) = -7t for all t > (72) 
v ' Jo t(?) ~ v ' 

The particle energy E = —Uip 1 = p° is 

[ R~ l (t) [{ Pi y + {p 2 f + (p 3 ) 2 ] 1/2 (m = 0) t < t 

E = (73) 

[ K + x-\t)( Pl f + Y-\t)(p 2 y + z-\t)( P3 y] 1/2 t>t 

with w 2 = (pi) 2 + (p 2 ) 2 + (p 3 ) 2 for all t > 0. (Note that w = R{t)E for t < t .) 
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The particle distribution function is given by 

f(t,Pu)=f(w) + h(p u )e- r ^ 

where 

- = f F(R(t)E) t < t 

1 \ F(t,E)+F I/ (t,E)e u + --- t>t 1 ' 

with F u ...{t ,E) = 0, F(t ,E) = F{R{t )E), and 

h( Pv )=H(t,E) + H v (t,E)e v + "- 

with H v ...(t,E) = H v ...(R(t)E)) for * < t . Note that in the FRW phase (f < t ) each 
F v ... and H v ... is a Liouville solution (hence functions of w only), because the multipoles 
decouple in L[f] = due to = 0. This is no longer true in the Bianchi I phase 
(t > t ) and hence F v ... and H v ... are functions of t and E. 

The distribution is specified to yield an Einstein-BGK solution for t > 0: 

H(t, E) = = H„(t, E) = F„(t, £) = = F VK (t, E) 

The resultant solution is 

/(*, p v ) = F(t, E) + e- r ^H Kp (t, E)e K e p + ■■■ (75) 
with F = F(R(t)E), H VK = H VK (R(t)E) for t < t , and H VK obeys 

H\(t, E) = diag [Vi (f, E), V 2 (t, E), -V x {t, E) - V 2 (t, E)\ 

where for t <t : 

V v (t,E) =V u {R{t)E) (76) 

/ w 3 VJw)dw = 
Jo 

and for t >t : 

Tt\(t) = ^e~ m / (E 2 - m 2 f 2 V u {t, E)dE (no sum) (77) 

15 Jm 

Finally, F and V v are subject for t > t to the constraint (67): 

- r (f-!)^(^)-il-!)^(^) 
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It is clearly possible to find F, V v such that (76) - (78) are satisfied (and therefore 
all field and conservation equations will be satisfied). The solution has zero number 
flux, bulk viscosity and energy flux for allt > 0. 

The early universe temperature is high enough that the particles have m — 0, and the 
matter distribution is described by (|75D, consistent with an FRW spacetime geometry. 
The conditions (67), (76) and (77) can be satisfied by appropriate H UK ^ and thus 
the field and conservation conditions are satisfied for both the FRW and Bianchi I 
phases. As the universe expands and cools below the threshold energy (at t = t ) for 
the particles under consideration, the distribution is no longer effectively of zero rest 
mass particles and the condition (76) is no longer satisfied. Condition (|28|) is required 
to ensure vanishing anisotropic stress in the FRW phase (i.e. the field equation -r vk = 
is satisfied). For the AW collision model (76) is satisfied by H UK ^ only if m — 0, by 
(p8|). Thus, as soon as the particle rest mass is no longer effectively zero, the anisotropic 
stress becomes non-zero (n UK ^ 0) and is given by (77). As a result shear anisotropy 
emerges, forcing the universe to evolve away from the FRW to the Bianchi I phase. 

During the phase t < to the collision rate is high and therefore e~ r may become 
small. This forces the non-equilibrium, anisotropic distribution function f(t,p u ) to 
approach the isotropic, collision-dominated equilibrium distribution f(w). However, as 
long as e~ r > (even if it is very small) the model presented here works. The model 
therefore represents a high temperature situation where the matter distribution is nearly 
isotropic (forced by the high collision rate). The distribution function has Bianchi I 
symmetry but satisfies all the conditions for an FRW universe. As the universe cools 
the particle threshold energy is reached and the FRW condition is no longer satisfied. 
The remnant of the initial anisotropy therefore acts as the seed for the change to Bianchi 
I geometry when the particle mass becomes significant, communicating the anisotropy 
of the distribution function to the spacetime geometry. Hence, this model suggests a 
mechanism by which anisotropy present in the radiation era could act as the seed for 
anisotropy generation in the spacetime geometry. The small anisotropy in the microwave 
background radiation could be a physical example of the remnant anisotropy considered 
here (see |16| for a related discussion). 
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